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\ Abstract 

We discuss the relation between matrix models and the Seiberg-Witten type (SW) theories, recently 
proposed by Dijkgraaf and Vafa. In particular, we prove that the partition function of the Hermitean one- 
matrix model in the planar (large N) limit coincides with the prepotential of the corresponding SW theory. 
This partition function is the logarithm of a Whitham r-function. The corresponding Whitham hierarchy 
is explicitly constructed. The double-point problem is solved. 

CD . 

1. It is well known that partition functions of matrix models are r-functions of integrable hierarchies of the 
Toda type [EJ. In the specific double scaling limit, these r-functions become r-functions of various reduction 
of the KP hierarchy If one makes the simplest, large- N (planar) limit, the partition function becomes 
the r-function of the dispersionless Toda hierarchy, which in turn becomes the r-function of the dispersionless 
(reductions of) KP hierarchy |j| after performing the continuum limit (which basically means working nearby 
a singularity of the partition function). All these dispersionless hierarchies are just Whitham equations over 
trivial solutions to integrable (Toda, KP) hierarchies. 

When solving matrix models, most attention was paid to one-cut solutions where the limiting eigenvalue 
distribution spans one interval on the real axis The results on multi-cut solutions || were few || |?| ||]. 
Recently, Dijkgraaf and Vafa proposed || the new insight on the multi-cut large- N limit of matrix models. 
Namely, they associated this limit with a Riemann surface and some related SW system. Its prepotential, 
which we prove here to be the logarithm of the large-A^ partition function, is typically associated with the 
logarithm of some Whitham r-functions [ll[] . This hints that the matrix matrix model in the large N limit 
of multi-cut type describes the Whitham system over a non-trivial, finite-gap solution to integrable (Toda, KP) 
hierarchy. In particular, this solution passes to a finite-gap solution of (reductions of) the KP hierarchy in the 
continuum limit. 

In this paper, we restrict ourselves with the simplest example of the Hermitean one-matrix model. We 
show that coefficients of the potential of the model gives rise to Whitham flows and manifestly construct this 
Whitham system. In fact, the authors of || associated the N = 1 SUSY gauge theory studied in with the 
SW system related to the multi-cut planar limit of matrix models. From the point of view of N = 1 SUSY 
theory, these coefficients must be identified with couplings in the tree superpotential, while the SW moduli 
are associated with v.e.v.'s of the gluino condensates. This gives an interpretation of the results of Jl2| in the 
Whitham hierarchy terms. 
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2. We call SW system the following set of data!: 



• a family Ai of Riemann surfaces (complex curves) C whose dimesnion coincides with the genus;u 

• a meromorphic differential dS whose variations w.r.t. moduli of curves are holomorphic. 

This data allows one to define the notion of prepotential |l3[ [l7| related to some integrable system |n| . 
Indeed, one can introduce variables 

cij = <j> dS (1) 

where Ai are A-cycles on C. Then, 

ddS 

dwi = — 2 

are canonical holomorphic differentials on C (normalized so that §. dujj — Sij). Then, introducing i?-cycles 
conjugated to A-cycles: Ai o Bj = Sij, where o means intersection, we obtain that 

d 

- I < ! - s = 

is the period matrix of C and is therefore symmetric. Hence, there exists a prepotential T such that 



— — d> dS = (p du)j = Tij (3) 
(JCLi ./ R . 



dS (4) 

oai J B . 

3. Let us consider the Hermitean one-matrix model. Its partition function is given by the integral over 
Hermitean N x N matrix M 

Z N = [ VMe- NttV W (5) 



where DM is the Haar measure on Hermitean matrices and the potential V(x) is a polynomial of degree n + 1. 
After integrating out angular variables, we obtain |Q 

In the large N limit, it is standard to introduce the density of eigenvalues 

i 

Then, (0) can be rewritten as 



Z 



N 



In the large N limit, this integral can be evaluated by the saddle point method. We then assume p(X) to be a 
continuous function such that 



p(X) > 0, / P (X)d\ = 1 (9) 

We then obtain the saddle point equation 

V{\) + £ = 2 / p(A') log(A - \')d\', A, A' e supp( /9 ) (10) 



where the support of the function p(X) comprises A such that p(X) ^ 0. It emerges in this equation because 
of the first condition in (||). The constant £ in equation ([To]) is just the Lagrange multiplier for the second 
condition in (|^). 

1 Various properties of such systems can be found in Q . 

2 Our definition of the SW prepotential does not imply any connection with prepotentials of N = 2 SUSY gauge theories 
Moreover, the prepotentials discussed in this paper are rather related to superpotentials of N = 1 SUSY theory [113] . 

3 This restriction can be waved, see examples in Jl6| . 



However, in order to use analytic tools (the Cauchy problem), we must investigate not eq. (|l0|) but its 
derivative 

V'(X) = 2 fj^dX', or V'(X) = / 4&jdX. (11) 



A- A' 



In order to solve this equation, we introduce the function 



(12) 



such that its imaginary part coincides with p(X) because of (|lT|). In the large N limit, it satisfies the equation 

V'(A)-V'(V) 



y'(X)-V"(X)+4 



A - A' 



'-p(X')dX' ee y 2 (A) - V' 2 (\) + / n _x(A) = 



where / n _i(A) is a polynomial of degree n — 1. This means that the general solution to ( p"T| ) is 

2n 

y 2 = V' 2 (X) - / n _x(A) = JJ(A - /ij) 

i=l 



(13) 



(14) 



This equation describes a hyperelliptic curve of genus n — 1. It is, however, not arbitrary for a fixed potential 
V(X), because it follows from ( [l3| ) that 



l/(A) - V"(A) = W(X) ~ - + O (A- 2 ) 



(15) 



A — *oo 



i.e., /ii are not independent. Here W(A) is the standard loop mean One more restriction comes from the 
normalization condition in (Q), and one is left with n — 1 free moduli. This is exactly what we need for SW 
system given on curve (|l4|). 

A solution to eq.([ll|) is parameterized by n — 1 moduli, i.e., these moduli span the moduli space of planar 
limits of the matrix model. The function p(A), which is imaginary part of y(A)/(27r), has the support on n 
different branching cuts. 

Note that within the standard matrix model framework, there are two more requirements that leave no 
moduli in solution. First of all, one can easily see that the sign of p(A) changes when coming to the next cut. 
This spoils non-negativity of p(X) and means that solution ( |l4| ) is not stable. This means that the cuts with 
p(X) negative must shrink to produce double points. Thus, y(X) becomes proportional, besides the square root 
of a polynomial, to some other polynomial that has odd numbers of zeroes between cuts and, therefore, changes 
sign on every next cut. 

We shall explain below that one can easily include these double points into the_.general SW and Whitham 
framework. Moreover, they allow one to construct more general Whitham systems.u 

The second requirement looks more fundamental. Namely, returning to original equation ([lO]), one has to 
check that the Lagrange multiplier £ is the same for every cut (while (ThO) only guarantees it is a constant on a 
cut). The difference of values of £ on two neighbour cuts is equal to H 



M2i + 1 



(16) 



where the integral runs from the right end of the left cut to the left end of the right cut. This gives n — 1 
additional constraints and leaves no moduli (there still remains a freedom in the number of cuts). We must 
wave this requirement in order to make the Whitham system nontrivial. So, instead of just matrix models, it 
is better to speak about matrix- model-like solutions of Cauchy problem (pi]). So, at the moment we just ignore 
this last restriction and work with solutions to eq.(|li"|), which we call matrix model solutions. We shall return 
to this point later. 



4. Now we associate an SW system with the planar limit of the matrix model. The family of genus g curves is 
described by eq.()l4|) with the restriction ([l5|). We describe it on the complex plane A by n = g + 1 cuts, Fig.l. 
Besides canonically conjugated A- and B-cycles, we also use the linear combination of -B-cycles: Bi = B^ — -Bj+i, 
B n -\ ee B n _\. Therefore, -B-cycles encircle the nearest ends of two neighbour cuts, while all B-cycles goes from 



4 V.Kazakov suggested to overcome non-stability of solutions without double points via some proper analityc continuation. 



a given right end of the cut to the last, n-th cut. For the sake of definiteness, we order all points /i, in accordance 
with their index so that jii is to the right of fij if i > j. 




Fig. 1. Structure of cuts and contours. 

The SW differential is 

dS = y(X)d\ (17) 
Its variations w.r.t. moduli is holomorphic on C ( |l3|) because all moduli are hidden in the polynomial /„_ i(A): 

ddS _ 0/ n _i(A) dX 



^moduli ^moduli y 



(18) 



This expression is holomorhic, because the leading coeffient of / n _i(A) is fixed by the normalization condition 
(^), and the differentials A fc ^ are holomorphic on the curve C ( |l3| ) for k = 0, 1, ...,n — 2. Therefore, we 
introduce the variables 

(H = -f y(X)dX = lm y(X)dX = p(X)dX, i = l,...,n-l (19) 

that have meaning of "the occupation numbers" (or numbers of eigenvalues) associated with a given cut. Note 
also that 

n-l 

; a t (20) 



\ I y(X)dX = Im f " y(X)dX = 1 - £ , 

1 J A„ Jn2n-1 i = i 



2 

which follows from (^|). This is exactly the condition that fixes the leading coefficient of /„_i(A) and leaves 
n — 1 moduli. It means that 

1 f 1 /* ddS 

- (p dwi = - (p — — = -1 for all i (21) 

2 JA n 2 JA„ » a J 

Now one defines the prepotcntial 

^ = f dS (22) 

OCLi J B . 

This prepotenatial is equal to logarithm of the matrix model partition function, logZ^ in.-the planar limit. 
Indeed, let T> be a set of contours A\ U U . . . U A g U A n . The large N partition function isa 

log Z N = - 1 £ y(X)V(X)dX + J/(A) log(A - A>(A')dA dV (23) 

We now calculate the derivative of logZ/v w.r.t. af. 

The expression in the brackets on the rhs of ( pi|) is a step function, which is equal to £j on each cut and its 
values on different cuts are (16) 



6 + is »(A')dA' = /ia for A e A 2 , 



£(A) = V(A) - / dA'log(A - A')y(A') = I 



£i = hi for A G Ax, 

: ( 2 5) 
& + i : B 1 us 2 u...us a _ 1 y( A ')^A' = ft s for A G A fl , 
& + /b 1 ub 3 u...ub s _ 1 ub s I/(A')dA' = &n for A G A„. 

5 From now on, we consider symbols <f> and res with additional factors (2iri)~ 1 so that reso^ = — resoo^ = <j> ^ = 1. 



We therefore have 



— ^p— — ■ = - / ^h(X) =-( ^(A) = -Zi + Zn =1 dS^l dS (26) 

and Zjy in the planar limit can be, indeed, identified with the prepotential e^. 



5. One can learn two lessons from this fact. First of all, we can return to the interpretation of different £i on 
different cuts within matrix model. The standard matrix model case of equal £j 's can be now formulated as the 
set of conditions 

f = ... = f=0 (27) 
oai oa g 

These are the conditions of minimum of the matrix model partition function w.r.L the occupation numbers. 
They can be removed by introducing different chemical potentials for different cuts.B We do not enter here any 
further details and go instead to another lesson. 

We know from studies of matrix models that their partition functions are r-functions of some integrable 
hierarchies JjJ. What are they in the planar limit? We have just proved that such a partition function is an SW 
prepotential in this limit. One typically associates logarithms of Whitham r-functions with SW prepotentials 
|To| , Oj. Therefore, we may expect that the matrix model partition function becomes the r- function of some 
Whitham hierarchy. An additional evidence for this comes from looking at the simplest one-cut large- TV solution 
of the matrix model, when the partition function becomes the r-function of the dispersionless Whitham hierarchy 
||. Now we construct this hierarchy in very manifest terms. 

First, we return to the problem of double points. Let us assume that some of the cuts shrink, i.e., 



y(X) =M„_ fe (A), 



2/r 



, ih) = M n - k (\)y/g 2 kW (28) 



where Af„_fc(A) is a polynomial of degree n — k and <72fc(A) is a polynomial of degree 2k. This means that one 
is effectively left with a new curve 

y(A) = Vg^W (29) 

This curve of lower genus k — 1 along with the differential dS = ilf„_fc(A)y(A)(iA remarkably give rise to a new 
SW system that depends on k — 1 moduli. 

To see this, one needs to take into account that there still holds eq. (|f3|) , 

y 2 (\)MZ_ k (\) = V'\\) - / n _x(A) ee V' 2 {\) 2(V'(X)W(X)) + , (30) 

where we let (•)+ denote the polynomial part of the expression in brackets. Then, varying dS and using j29|), 
we obtain for the general variation SdS: 

SdS = 5 (M n _ fc (A M A)) dX = WA)*M„- fc (A) + lM„- fc (A)^A) dA 

y 

On the other hand, doing a variation 5 of M n _k(X)y(X) that does not alter the potential, we obtain from ( |30| ) 
that 

~ SdS = ~\ J fn m X L dX - ^ 

2 M n _fc(A)y(A) 

Because this variation is a particular case of (|3l|), we obtain that zeroes of M n ^k(X) in the denominator of 
( |32| ) must cancel, so the maximum degree of the polynomial in the numerator is n — 2. The variation is then 
holomorphic on curve (p9|). 

This solves the problem of double points. The corresponding system with double points (the large- TV limit 
of the matrix model) is still described by the SW theory. 



6 Putting differently, one can interpret these conditions as a criterium of stability against tunneling of eigenvalues between 
different cuts Stability is achieved by imposing equality of the chemical potentials of all cuts. 



6. Let us return to the case n = k. This SW system is described by n — 1 and n + 1 additional parameters. 
They are expected to be Whitham times giving flows on the moduli of the finite gap solution. In this moduli 
space there are deformations leaving the curve within the family (variations of /„_i(A) that does not change 
the potential), and those transversal to the family. These latter are defined by the potential. If one stays within 
the family of hyperelliptic curves (p"3|), there are exactly n + 1 transversal deformations. Therefore, in order to 
have enough many deformations (Whitham times), one needs to involve potentials of high enough degree, i.e. 
to deal with the construction with double points. 

In our manifest construction of the Whitham system we mainly follow Jll], (see also ^l)). In order to 
construct a Whitham system, one needs to add to the SW data a set of punctures with local coordinates in 
their vicinity. These points here are the two infinities on the curve ( p9[ ) and the local parameter is rj — j. Now 
one introduces a set of meromorphic differentials dQ n with the poles only at punctures (since the hyperelliptic 
cur ve (|29|) is invariant w.r.t. the involution y — > —y, from now on we just work with either of the two infinities, 
see |20{) and the behaviour 

dQ m = {r]-" 1 - 1 + 0(1)) dq, f] -» (33) 

Then, the Whitham system is generated by a set of equations for these differentials and the holomorphic 
differentials dut;: 



ddflp ddQ. m ddVL m ddwi dduji ddiOj 



dt m dtp ' dcii dt m ' daj dai 



(34) 



These equations implies that there exists a differential dS such that 



dds dds 

— — = dwi, — — = dilm (35) 

OCLi ot m 

Let us check that the differential dS = M n _k(X)y(X)dX given on the curve ( p9| ) with the relation for moduli 
@ really satisfies @. 

Indeed, we have proved the first set of relations ( p5| ) in the previous paragraph. Now let us consider variations 
of the potential, i.e., variations w.r.t. Whitham times t m . Then, we obtain instead of (132 



2 M n - k (X)y(X) V ' 

while (|3l|) still holds. Repeating the argument of the previous paragraph, we conclude that the zeroes of 
M„_fe(A) cancel from the denominator and, therefore, the variation may have pole only at A = oo or r/ = 



0, i.e. at the puncture. In order to estimate this pole, one needs to use (31), which implies that dS = 
M n _k{X)y(\)d\ — > (V'(X) + 0(j))dX and, therefore, the variation of dS at large A is completely determined 
by the variation of V'(X). Parameterizing V(X) — one comes to ( |35| ) up to a linear combination of 

holomorphic differentials. One may fix the normalization of dQ m that are also defined up to a linear combination 
of holomorphic differentials so that eq. ( |35| ) would be exact. What does this normalization mean? Throughout 
all our consideration we deal with a^s and t^s as independent variables. This unambiguously defines the way 
a^'s depend on the coefficients of and is achieved merely by imposing the (obvious) condition Jll], ||o| 

da.. 



7 dfL m = 0Vi,m, (37) 

Thus, similarly to eq.(Q) we can invariantly introduce variables t m via the relation 

t m = res n=0 i] m dS (38) 
and define the prepotential that depends on both at and t m via the old relation (Q) and the similar relation 

|^ = -rc &n=0 rr m dS (39) 
ot m m 

One can immediately prove that such a prepotential exists Jll], |20), i.e., the second derivatives are symmetric, 
and, moreover, similarly to Sec. 4, we find that thus defined T coincides with logZ/v in the planar limit. For 
this, we apply the formula similar to (|24|) with the only difference that the potential V^A) itself is changed. We 
then obtain from ( p5| ) 

^ ~5 ■ ^ - \ L d ^ - - i - " J + i^-^ w 

which by virtue of ([37]) gives ([39|). Therefore, Zn in the planar limit is the Whitham r- function, and the whole 
machinery of Whitham systems works here in full strength. 



7. After having constructed the large N (planar) limit, next step is to take the (double scaling) continuum 
limit. Namely, one has to dwell nearby a singularity (branching point) of p(A). Say, one can work nearby the 
left end of the very right, n-th cut /i2n-l (§]■ The standard argument then is that one feels no other cuts, since 
they are far away. This would mean that multi-cut solutions coincide with the one cut solution in the continuum 
limit. This is, however, the case only if all the other branching points do not come close to /U2n-i- Otherwise, 
there exist non-trivial continuum limits Q|. In the most non-trivial situation, all the branching points but H2n 
come close to each otheiu. This is equivalent just to sending X 2n to infinity. Such a curve describes a finite 
gap solution to the KdV hierarchy, moreover, the corresponding SW system is also associated with KdV p2| . 
Therefore, we expect that the matrix model partition function in this limit describes (in leading order) the 
Whitham hierarchy over KdV finite gap solution. 

However, it would be very instructive to construct an entire double scaling limit in this situation, in partic- 
ular, to fix proper scaling behaviours. This means to match properly the growth of N and approaching to the 
singularity. Then, one could address the problem of exact (matrix model) <-> (SW) correspondence, in particular, 
in integrable terms. In particular, it would be interesting to see what is the proper deformation of SW systems 
in this case. 

In this respect, the problem of studying higher-genus corrections looks very natural because the whole 
matrix- model-like solution (the solution to the loop equation, see 0j) in all genera is completely determined 
by the set of data {cii,t m } and because it was proved [p3j that spectral correlators manifest the universality 
property for multi-cut solutions as well. The integrable system that will appear in this approach must be a 
generalization of a Whitham system. 

At last, let us note that the construction considered in this paper is directly extendable to other matrix 
models. In particular, one can consider the model of normal matrix that has much to do with the problem of 
Laplacian growth pj| . In fact, its naive large N limit describes how external and internal moments of a domain 
are related. This domain is a counterpart of the one cut. Moreover, the system is described by the Whitham 
hierarchy that is the dispersionless Toda system. Now, considering a multi-domain solution and introducing 
chemical potentials for different domains, one has to get the Whitham system over a finite-gap solution to the 
Toda system. This Whitham hierarchy should relate the external and internal moments of several domains. 

The authors are grateful to A.Gorsky, S.Gukov, V.Kazakov, A.Marshakov and especially to A.Morozov for 
fruitful discussions. We are also due to Chung-I Tan for providing us with proper references on the multi-cut 
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